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Abstract. The 2nd order corrections are obtained to both forward and nonforward interaction of 
reggeized gluons in the octet colour channel using as a basis the bootstrap relation and a specific 
ansatz to solve it. The obtained forward kernel coincides with the logarithmic term plus two first 
non-logarithmic terms in the Pomeron 2nd order kernel. Both forward and nonforward kernels are 
found to be infrared finite. 
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1 Introduction. 

Recently the 2nd order corrections were calculated for the forward hard (BFKL) pomeron equation 
[1,2]. They clearly divide into a correction for the running of the coupling and the 2nd order corrections 
in the running coupling proper. Generalization of this result for the other physically important case, 
that of the odderon, requires knowledge of the 2nd order corrections to the non-forward interaction 
of two reggeized gluons in the octet colour state. It is well-known that this interaction is severely 
restricted by the so-called bootstrap relation, which is, in fact, the unitarity requirement for the octet 
channel [3]. The form this relation takes in the 2nd order was recently discussed in [4]. 

In previous publications one of the authors (MAB) proposed to use the bootstrap relation and 
a specific ansatz to satisfy it for the introduction of the running constant, to all orders of the fixed 
one, into the gluon octet interaction [5]. Since in the lowest order this interaction is just one half of 
the interaction in the vacuum (pomeron) channel, this method also serves to find corrections for the 
running coupling in the pomeron equation both for forward and non-forward cases. Note that whereas 
for the forward case these corrections are trivial and can easily be recostructed via the renormalization 
group, this is not so for the non-forward case, with more than one scale. 

In this paper, using the bootstrap and the ansatz of [5], we calculate the full second order correc- 
tions for the non-forward gluon interaction in the octet channel necessary for the odderon equation. In 
this way we also automatically find the corrections for the running of the coupling for the non-forward 
pomeron equation. 

2 Basic equations 

Up to the 2nd order in the fixed coupling, two gluons are described by a Schrodingcr equation in the 
transverse space 



Here q\ and q 2 are the momenta of the two gluons; E = 1 — j where j is the angular momentum. 
The Hamiltonian H is a sum of the kinetic energy, given by a sum of gluonic Regge trajectories with 
a minus sign —u{q\) — Lo{q 2 ), and the interaction — V^ R \ which depends on the colour state R of the 
two gluons. We shall be interested in either the vacuum state (R = 1) or the octet state (R — 8). As 
mentioned, to the lowest order in the running coupling, V^ 1 ) = 2V^ . 

The bootstrap relation requires that the solution to Eq. (1) for the octet channel be the reggeized 
gluon itself. This can be fulfilled if the interaction V$ is related to the trajectory as 



H<j)(quq 2 ) = E(f>(q 1 ,q 2 ). 



(1) 





(3) 



and 




viv'i) v(q' 2 ) v{qi-q'i) v(q[)v(q' 2 ) 



v(qi) , v{q2), 1 v(q) 



) 



(4) 
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In the lowest order of the fixed coupling one has (we take N c = 3) 

The forms (3) and (4) guarantee that to the lowest order in the running coupling the full interaction 
in the vacuum channel is infrared stable [5]. A stronger result is that they also guarantee that the 
sum of single and pair terms in the odderon equation stays infrared stable in all orders in the running 
coupling. Indeed, this sum can be represented as a sum of three pair terms — (1/2) (K\2 + K 2 3 + if 31 ) , 
where e.g. 

K 12 = w(q 1 )+w{g 2 ) + 2V^\ (6) 

and this combination is infrared stable provided 77(g) goes to zero as q 2 — > modulo logarithms. Of 
course, this does not settle the question of the eventual infrared stability of the odderon equation in 
the higher orders, since already at the 2nd order triple interaction terms seem to appear. We shall 
not deal with this problem here. 

Our aim will be to study the 2nd order corrections to the kernel K which follow from the 2nd 
order corrections in 77(g). The latter can be found from the known form of the trajectory u)(q) to the 
2nd order in the fixed coupling [6]: 

uj(q)=J 1 \q)+w^(q), (7) 

where the superscripts refer to the order in the coupling, and, in the dimensional rcgularization in the 
MS scheme, 

u>W{q) = -f + 2t + e(t 2 - 27/(1))) , (8) 

^ (2) (q) = -t {m± - t 2 ) + B(- e + 2t) + Cj , (9) 

Here and in the following t = ln(g 2 //i 2 ); \x is the normalization point; e — > 0; 

_ 2 _ 3a s r(l - e) 
9 - ( 47r )i+ £ 

and g 4 is just the square of (10). The coefficients A, B and C are 



(10) 



We have retained a term proportional to e in (8) for future reference. 



3 Function 77(g) to the 2nd order 

We present function 77(g) in the form 

77(g) =r/°)(g)(l+£(g) + ...), (12) 

where £ represents the terms of the 2nd order in a s . The substitution of the running coupling instead 
of the fixed one into (5) evidently gives a contribution to £ 

T^- 
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This term is what we call a correction for the running of the coupling. However there may be other 
terms in £ which represent corrections of the 2nd order in the already running coupling constant. We 
are going to find all these by matching the form of lu^ found via function r)(q) given by (12) with the 
expression (9) found from the perturbative QCD. 

Putting (12) into Eq. (3) we obtain in the 2nd order 

- ( %)=" ( %K(?) + ^/^«- (13) 

With uj^(q) known, this is an integral equation for £. We shall solve it by choosing an appropriate 
functional form of £(</): 

i(q) = c + dt + eft\ (14) 

where the coefficients c, d and / in principle may depend on e but are regular at e = 0. We introduce 
an extra factor e into the third term of (14) to avoid poles of the third order in e absent both in lj^ 
and the first term on the r.h.s. of (13). 

We denote the integral term in (13) as X. In the dimensional regularization we have 

X = ^t(cX + dX 1+ feX 2 ), (15) 



= J^W> (16) 



where 

X 

D = 2 + 2e and t\ = \n(qf/^ 2 ). X n can be calculated as the n-th derivative in a, taken at a = of a 
general integral 

d_Jiq_ , , 

This integral can be easily calculated by standard methods using the generalized Feynman represen- 
tation. We find: 

1 = 77 {q) T(l-a)T(a + 2e) ' (18) 

^From this we obtain 

Xo=* (q) ^ , (19) 

Xi=X a Y, y = (t-V(l-e)+V(e)+V'(l)-V'(2e)), (20) 

X 2 = X Q (Y 2 + 1/(1 - e) + V'(e) - - ^(2e)). (21) 



One has to take into account that the 1st order trajectory (8) is related to X as 

- ( %) = -tS4*o. (22) 



Combining the two terms on the r.h.s of (13) we get 



" (2 H q ) = j^Xo(-- e +c + ht + fet 2 ) (23) 

with u = d — 2/. Recalling the form of u>^\ Eq. (8), separating g 2 and neglecting terms of the order 
e or higher, we finally get 

J 2 \q) = 2f (-« + ^ +ci + \dt 2 • (24) 
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Comparing this expression with the one found from the perturbative QCD (9), we note at once 
that to match them we have to assume that both u and c contain terms of the higher order in e: 

u = u Q + me + u 2 e 2 , c = c + C\t. (25) 

With this form the 2nd order trajectory becomes 

^ 2 \q) = 2.9 2 (-^ + ^-^ + cot + \dt 2 + Cl - u 2 + Uo V'(l)) (26) 

. So we have just 5 parameters Co, Ci, d, Uo, Ui, u 2 to match this expression with 5 terms in (9) 
containing 1/e 2 , l/e,t, t 2 and a constant. We get for our parameters 

w = (1/2) >f : A, co - «i = -(l/2)g 2 B, c - -g 2 B, d = g 2 A, u tp'(l) + a - u 2 = -(1/2) g 2 C. (27) 

We observe that to match the perturbative trajectory it is not sufficient to take into account only the 
corrections to 77(g) due to the running of the coupling, which corresponds to all coefficients equal to 
zero, except for equal d and u . Quite a few proper 2nd order corrections have to be also included. 



4 Forward interaction in the octet channel 

Having determined the 2nd order corrections for 77(g) we are now in a position to find these corrections 
for the interaction of gluons in the octet channel. As mentioned, to the lowest order in the running 
coupling, the vacuum channel interaction is just twice the octet one. So we also find the corrections 
due to the running of the coupling for the vacuum channel. In this section we calculate the 2nd 
order interaction for the forward case q = 0, which is a much simpler task due to the fact that the 
eigenfunctions of the 1st order kernel are simple. 

As for the pomeron, at q — the eigenvalue equation for the octet channel interaction 

— K(f> = E<f>, (28) 

with K given by (6), can be simplified by substituting 

(j){q) -> T}(q)<j>(q). 

Then (28) takes the form 

+ £sf - - W (29) 

In the 1st order this is just the BFKL equation. In the 2nd order we get for the left-hand side 

4tt 2 J (q- q'y 

where £ is defined by (12) and determined by (14) and (27). 

We shall calculate the action of the kernel on the 1st order BFKL eigenfunctions 

MQ) = V 2X - (31) 

In fact, for the normalizability one should take A = — 1/2 + iv and multiply (31) by the appropriate 
normalization factor, but this is irrelevant for the properties of the kernel K^ 2 \ 
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We again introduce a notation X for the integral part of the left-hand side (without the minus 
sign). Putting (14) into the integrand we find the same form (15) for X, with an extra factor two and 
X n which are now given by 

-r|H|< 92 = Q - qi, h = Mqi/H 2 )- (32) 
Qi qi 

These integrals can again be obtained by differentiating in a at a = the generic integral 

f d D q iq 2a 1+ 2 +x T(-\-a-e)T(a + e)T(l + \ + e) 
J q^ 2X q 2 {q ' T(l-a)T(-\)T(l+a + +X + 2e)- 

^From this we find 

y ^i+ £rn 2^ + A r(-A- £ )r( £ )r(i + A + e) 

Xa = n {q) r(-A)r(i + A + 2 £ ) ' (34) 

X 1= X Y, F = i-V(-A-e)+^(e)+V(l)-V'(l + A + 2 e ), (35) 
X 2 = X a {Y 2 + <//(-A - e) + V'(e) - /(I) - V>'(1 + A + 2e)). (36) 
Combining the terms we get for X 

- +c + ?i(i + A)+ 

e(d(V'(-A) + - 2^(1 + A)) + /((t + A) 2 - ^'(-A) - 2/(1) + 3^(1 + A)))] , (37) 
where we denoted 

A = 2^(1) - - V(l + A) (38) 

the eigenvalue of the 1st order kernel (up to a constant factor). To obtain our final result we separate 
a factor to be included into g 2 and develop the rest, including X , in e up to constant terms. Then 
we get 



v _ 12t* s 



X = 4gy A J + °- + c(t + A) + ld((t + A) 2 + V'(-A) - ^'(1 + A)) + u^(l)) . 



(39) 



Now we have only to subtract it from the 1st term on the righthand side of (30). Taking cuW in 
the form (24) we find that the terms singular in e cancel, so that the kernel turns out to be infrared 
stable, as expected. Our final result for it is 

KPty = Ag 2 (-dAt - cA + ^d(-A 2 + ?//(l + A) - V'(-A))) ip(q). (40) 

One notes that out of five coefficients entering £ only two remain in the final expression for the action 
of the kernel. Using (27) we find 

K( 2 V = 4(^) 2 {-AAt + BA + ^A(-A 2 + V>'(1 + A) - V'(-A))) ^(q). (41) 

The first term, proportional to \n(q 2 / fj, 2 ) gives the correct contribution from the running of the 
coupling. As mentioned this term is the same for the octet and vacuum channels. So the bootstrap 
indeed gives the correct description of the running of the coupling, at least, in the 2nd order. Other 
terms are corrections to the interaction in the 2nd order in the running coupling constant. They 
refer only to the octet interaction since at this order the vacuum channel interaction is different. It 
is remarkable however that they exactly coincide with the first two terms found for corresponding 
correction in the vacuum channel. As a whole the found expression is much simpler than the one for 
the vacuum channel, which includes many more terms with a complicated dependence on A. 



7 



5 The non-forward pair kernel in the octet colour state 

The bootstrap gives a relation between the gluon trajectory and its pair interaction in the octet 
colour state. As mentioned this interaction is important for the the odderon, where it enters in the 
combination (6). In this section we shall obtain the 2nd order correction to the kernel (6) using the 
found form of the function rj and Eqs (3) and (4) . 

It is convenient to present the kernel (6) in the form which explicitly shows its infrared stability. 
We use an identity 

1 1 1 



V(n)v(l2) v(Qi)(v(Qi) +77(172)) t)(Q2)(v(Qi) + V(Q2)) 
to present the trajectory (3) in an equivalent form 



(42) 



<"(«) = 2 I ,0 v r^TTTZ i VV 1 = * + 92- (43) 
J (2tt)- 2 ?7(< ? i)(j7(< ?1 ) +77(92)) 

Then denoting qn> = qi — q[ for brevity and combining the terms with r](qi — q[) in the denominator 
in the trajectories and in the interaction, we obtain for the action of the kernel (6) on an arbitrary 
function ip(qi) 

h r= , ; <*Vi»7(«i) ( Wi) V>(9: 



(27r) 2 r7( gi - q[) V V(q[) v(q[) + v(qi - «i) , 

+ a „ 2) - 2 [ fM^Ml (u) 

As in the case of the standard BFKL equation, this form shows explicitly that all singularities at 
small momenta q[, q' 2 and qw cancel and the kernel is infrared stable if function 77(g) goes to zero not 
essentially faster than q 2 . 

Now, preserving the infrared stability, we present in (44) function 77 in the form (12) keeping terms 
up to the 2nd order in a s . For the 2nd order contribution we obtain 

^ (2 V = 6« s / /f^(^(£(?i) ~ ^ ««iiO)- 
J y 2lT ) Qw q[ 

Qi +qfi> Qi +Qw Qi +Qw 

+ (1 - 2) - 6a s f f/\ q X iq :l m - (45) 

J (2irYq\ qL 



Presenting 



a' 2 a 2 

^ -l-^T^ (46) 



q[ 2 + q\ v Qi + in 



we separate from the first term a contribution 

j2 „i „2 



6a *) / TT^TT^-^Qi) - £{qu>)), (47) 
J (2TTy(q[ + qf 1 ,) 2 

which is equal to zero, since it changes sign under the substitution q[ ^ qi — q[. 
This brings us to a simple expression 

K m * - 6c " / ?Sf - +?r+) <««> - - ««■*» 

J WQw \ q[ q[ + q{ v ) 
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+ (1 - 2) - 6a s [ ^X^l m - (48) 

Now we have only to put our expression (14) for £(q) with the coefficients given by (25) and (27). 
Since the kernel is infrared finite, all terms in (14) proportional to e or e 2 give no contribution, 
and, as in the forward case, the result only depends on c = —g 2 B and d — g 2 A. It is also clear that 
the constant term Co in just rescales the first order kernel So in the end we finally obtain 

an explicit form for the action of K ^ as 

= ^1 \ BK (^ + 6Aas ( ( _ -4^lL-) In ■ ^ 



+ (1 - 2) - 6Aa s [ ^ In 4^1 • (49) 
J (2nyq[ 2 q> 2 2 q' 1 2 q' 2 \ 

As in the forward case, one expects that acting on the eigcnfunctions of the lowest order kernel, 

this 2nd order kernel will give a shift for the eigenvalue plus some terms which come from the running 

of the coupling. As mentioned, the latter should be the same as for the non-forward kernel in the 

vacuum colour state. Unfortunately, although the found kernel does not look too complicated, we 

have not able to find its action on the first order eigenfunctions up to now. Work in this direction is 

in progress. 



6 Conclusions 

We have found the 2nd order corrections to the interaction of two reggeized gluons in the octet 
channel, using the bootstrap relation and a specific ansatz to solve it proposed in [5] . This interaction 
is important for the odderon. The relevant kernel for a pair of gluons is found to be infrared finite 
and have a rather simple form. For the forward case, apart from an expected term which describes 
the running of the coupling, it contains a correction which coincides with the first two terms found 
for the vacuum channel in [1,2]. 

The infrared stability of the found octet pair kernel, if correct, implies that the three-body in- 
teraction in the odderon has to be infrared finite by itself. In our approach this stability, in fact, 
automatically follows from the applied ansatz to satisfy the bootstrap relation. We have no proof 
that this ansatz is unique. Calculation of the three-body interaction for the odderon and study of its 
infrared properties is thus essential to support our approach. 

After this work had been completed and published as a preprint [7] a paper by V.Fadin et al. [8] 
appeared in which the quark contribution to the interaction in the octet colour channel was directly 
calculated from the asymptotics of the relevant Feynman diagrams. It was claimed in [8] that their 
results disgagreed with ours and that therefore our bootstrap condition and the ansatz to solve it were 
wrong. In fact this claim is absolutely unfounded: the quark contribution to the interaction found 
by our method identically coincides with the results of [8] (see Appendix), although the amount of 
labour necessary to obtain it is incomparably smaller. 
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8 Appendix 

The authors of [8], in all probability, have not properly taken into account that the unsymmetric 
kernel V calculated in this paper is a different quantity as compared to their symmetric kernel K r . 
It is however straightforward to relate them. First one goes over from V to a symmetric kernel W in 
the trivial metric (that is with the integration measure d 2 q\) by expressing 



On the other hand, the kernel K r of [8] is defined in the metric with an extra factor q\q\ in the 
denominator. Taking this into account we obtain the desired relation 



Now we are going to show that, first, the quark contribution to the gluon Regge trajectory has the 
form implied by our anzatz (3). We shall determine the part of r/(q) coming from this contribution. 
Then using (4) and (51) we shall find the corresponding part of the irreducible kernel K r and compare 
it with the results of [3] . To simplify the comparison with [8] we shall use the relevant expressions in 
the unrenormalized form and for an arbitrary number of colours N . 

The part of the gluon trajectory which comes from quark is given by the expression [8] 




(50) 




(51) 




(52) 



where 



b 



g 4 NN F r(l-e)T 2 (2 + e) 
(27r) B - 1 (47r) 2 + e er(4 + 2e)' 



(53) 



On the other hand using (3) and (12) 




(54) 



where 



g 2 N 



(55) 




= -— <l' 



(56) 



Now we pass to the irreducible kernel K r . From (4) and (51) we express it via £(q) as 




V 9i 12 
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i\/T^«(«») + €(^)-«9i)-^)-2««nO)- 
V 92 <7i 



2 2/2/2 
9l<?29l <?2 



(2^)-€(«i)-€((ft)-^)-^)) 

(57) 



Putting here the found Eq. (56), we obtain the quark contribution 



1ll2 ,„2e , „/2£ „2e „/2e 
k 2 



I 2e 1 /at / rw . 



4C(9l e + «i* " " 92 2£ " 2* 2e ) - q 2 (2 q 2 < - qf - qf - q^ - qf) 



(58) 



Comparing this expression with the one found in [8] (Eq. (47) of that paper) we observe that they 
are identical. This means that our bootstrap condition and the anzatz to solve it are valid at least for 
the quark contribution in the next-to-leading order. 
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